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We present a search for a neutral, long-lived particle L that is produced in e + e“ collisions and 
decays at a significant distance from the e + e _ interaction point into various flavor combinations 
of two oppositely charged tracks. The analysis uses an e + e _ data sample with a luminosity of 
489.1 fb -1 collected by the BABAR detector at the T(4S), T(3S), and T(2S ) resonances and just 
below the T(4S). Fitting the two-track mass distribution in search of a signal peak, we do not 
observe a significant signal, and set 90% confidence level upper limits on the product of the L 
production cross section, branching fraction, and reconstruction efficiency for six possible two-body 
L decay modes as a function of the L mass. The efficiency is given for each final state as a function 
of the mass, lifetime, and transverse momentum of the candidate, allowing application of the upper 
limits to any production model. In addition, upper limits are provided on the branching fraction 
B(B —> X S L), where X s is a strange hadronic system. 

PACS numbers: 13.66.Hk, 14.80.-j, 14.80.Ec 


Recent anomalous astropliysical observations flf- 
liave generated interest in GeV-scale hidden-sector states 
that may be long-lived H£l2- Searches for long-lived 
particles have been performed both in the sub-GeV m3- 
[l5j and multi-GeV 1(1- 211 mass ranges. Dedicated ex¬ 
periments to search for such particles have been pro¬ 
posed [111 or are under construction 231. However, 
the 0(1 GeV/c 2 ) mass range has remained mostly unex¬ 
plored, especially in a heavy-flavor environment. B fac¬ 
tories offer an ideal laboratory to probe this regime. Pre¬ 
viously, the only R-factory results were from a search for 
a heavy neutralino by the Belle collaboration [24]. 


We search herein for a neutral, long-lived particle L , 
which decays into any of the final states / = e + e _ , R + , 
7r + 7r - , K + K~, or K±ttY . A displaced vertex and 
two-body decay kinematics constitute the main means 
for background suppression, and the search is performed 
by fitting the distribution of the L-candidate mass. 


The results are presented in two ways. In the “model- 
independent” presentation, no assumption is made re¬ 
garding the production mechanism of the L. Rather, we 
present limits on the product of the inclusive produc¬ 
tion cross section a(e + e~ -A LX), branching fraction 
B{L -A /), and efficiency e(/) for each of the two-body 
final states /, where A' is any set of particles. As supple¬ 
mental material to this letter [25|, we provide tables of 
the efficiency as a function of L mass m, transverse [36] 


momentum pp in the center-of-mass (CM) frame, and 
proper decay distance cr, assuming the L to be a spin- 
zero particle. The provided upper limits, efficiencies, and 
Pt distributions of the simulated events used to obtain 
the efficiencies facilitate the application of the model- 
independent presentation of the results to any specific 
model of L production. In the “model-dependent” pre¬ 
sentation, we provide limits on the branching fraction 
for the decay B -A X S L, where X s is a hadronic sys¬ 
tem with strangeness — 1. This presentation is motivated 
by Higgs-portal models of dark matter and other hidden 
sectors (8- 111. 


The data were collected with the BABAR detector at 
the PEP-II asymmetric-energy e + e _ collider at SLAC 
National Accelerator Laboratory. The sample consists 
of 404.0 ± 1.7 fb~ 1 collected at a CM energy correspond¬ 
ing to the T(4S) resonance, an “off-resonance” sample 
of 43.74 ± 0.20 fb _1 collected about 40MeV below the 
T(45) peak, 27.85 ±0.18 fb” 1 collected at the T(3S), and 
13.56±0.09fb _1 taken at the Y(2S) [26[. The T(4S) sam¬ 
ple contains (448.4 ±2.2) x 10 6 BB pairs, and the Y(3S) 
and T(2S) samples have (121.3 ± 1.2) x 10 6 Y(3S) and 
(98.3 ± 0.9) x 10 6 Y(2S) mesons, respectively [27]- An 
additional Y(4S) sample of 20.37 ± 0.09 fb -1 is used to 
validate the analysis procedure and is not included in the 
final analysis. 


The BABAR detector and its operation are described 
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in detail in Refs. [28| and [29j. Charged-particle mo¬ 
menta are measured in a tracking system consisting of 
a five-layer, double-sided silicon vertex detector (SVT) 
and a 40-layer drift chamber (DCH), both located in a 
1.5 T axial magnetic field. Electron and photon energies 
are measured in a CsI(Tl) electromagnetic calorimeter 
(EMC) inside the magnet coil. Charged-particle identifi¬ 
cation (PID) is performed using an internally reflecting, 
ring-imaging Cherenkov detector, as well as the energy 
loss measured by the SVT, DCH and EMC. Muons are 
identified mainly with the instrumented magnetic-flux re¬ 
turn. 

Using Monte Carlo (MC) simulations, we determine 
both the signal mass resolution and reconstruction effi¬ 
ciency. The events are produced with the EvtGen [3ll ] 
event generator, taking the A spin to be zero. We gen¬ 
erate two types of signal MC samples. In the first type, 
which is used to create the efficiency tables [25[ for the 
model-independent presentation, the A is produced at 
11 different masses, m^ c = 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 
9, and 9.5GeV/c 2 . For nig 10 < 4GeV/c 2 , the A is cre¬ 
ated in the process e + e _ — > BB, with one B meson de¬ 
caying to A + Nn ( N = 1, 2, or 3) and the other B 
decaying generically. At higher masses, the production 
process is T(4S') —>• A + Nn. In both cases, the A is pro¬ 
duced uniformly throughout the available phase space, 
with an average transverse decay distance of 20 cm. The 
events are subsequently reweighted to obtain efficiencies 
for other decay lengths. Note that these specific processes 
do not reflect preferred hypotheses about the production 
mechanism, nor do the results depend on these processes. 
Rather, they are a convenient method to populate the 
kinematic range for the efficiency tables. 

The second type of signal MC, used for the model- 
dependent presentation of the results, contains B —> X S L 
decays, for the seven mass values m = 0.5, 1, 2, 3, 3.5, 
4, and 4.5GeV/c 2 . The X s is nominally taken to be 10% 
K , 25% AT*(892), and 65% A'*(1680) j3C|, with the high- 
mass tail of the X s spectrum suppressed by phase-space 
limitations, especially for heavy A states. This choice of 
X s composition results in an A-momentum spectrum as 
a function of m* that reproduces the dimuon spectrum 
for B —> X s fi + p,~ in events generated with EvtGen using 
the BTOXSLL model 311. The other B meson in the 
event decays generically. 

In addition to the signal MC samples, background 
MC samples are used for optimizing the event selec¬ 
tion criteria and studying the signal extraction method. 
The background samples are e + e _ —>• BB (produced 
with EvtGen 31|), t + t~, pA (KK2f [32}), e + e - 
(BHWIDE [13), and qq events (JETSET Q) , where q is 
a u, d , s, or c quark. The detector response is simulated 
with GEANT4 Eal. 


The A candidates are reconstructed from pairs of op¬ 
positely charged tracks, identified as either e + e _ , p + p~, 
e^pN, 7r + 7r _ , K + K ~, or A' ± 7 t^. The PID efficiency 


depends on the track momentum, and is in the range 
0.96 — 0.99 for electrons, 0.60 — 0.88 for muons, and 
0.90 — 0.98 for kaons and pions. The pion misidentifi- 
cation probability is less than 0.01 for the electron PID 
criteria, less than 0.03 for the muon criteria, and aver¬ 
ages at 0.06 for the kaon criteria. A track may have 
different PID assignments and may appear in multiple 
pairs. Each track must satisfy do/ay 0 > 3, where do is 
the transverse distance of closest approach of the track 
to the e + e _ interaction point (IP), and oy 0 is the do un¬ 
certainty, calculated from the SVT and DCH hit position 
uncertainties during the track reconstruction. The two 
tracks are fit to a common vertex, and the y 2 value of 
the fit is required to be smaller than 10 for one degree of 
freedom. The two-dimensional vector f between the IP 
and the vertex in the transverse plane must have length 
r = |r| in the range 1 < r < 50 cm, and the uncer¬ 
tainty on r is required to satisfy ay < 0.2 cm. We require 
the angle a between f and the A-candidate transverse- 
momentum vector to satisfy a < 0.01 rad. The uncer¬ 
tainty a m on the measured A-candidate mass to must 
be less than 0.2GeV/c 2 . The A candidate is discarded 
if either of the tracks has SVT or DCH hits located be¬ 
tween the IP and the vertex, or if the vertex is within 
the material of the beampipe wall, the DCH support 
tube, or the DCH inner cylinder. Candidates must sat¬ 
isfy the following decay-mode-specific invariant-mass cri¬ 
teria: TO e + e — > 0.44GeV/c 2 , m /J + /J - < 0.37GeV/c 2 or 
to^+j^- > 0.5GeV/c 2 , TOg± M T > 0.48GeV/c 2 , m n + n - > 
0.86GeV/c 2 , m K + K - > 1.35GeV/c 2 , and m K ± n T > 
1.05GeV/c 2 . These criteria reject background from 
AT° —>• 7r + 7r - and A —> pn~ decays. In addition, other 
than in the p + p~ mode, they exclude low-mass regions 
in which the mass distributions of background MC events 
are not smooth, and therefore incompatible with the 
background description method outlined below. We re¬ 
quire at least one of the tracks of A —>• p + p~ candidates 
with to > 8GeV/c 2 to have an SVT hit. This rejects 
candidates that decay into p + p~ within the material of 
the final-focusing magnets, and thus have poor mass res¬ 
olution. These selection criteria are found to yield near- 
optimal signal sensitivity given the broad range of m and 
r values of this search. 

For each decay mode, we determine the full efficiency 
e, including the impact of detector acceptance, trigger, 
reconstruction, and selection criteria, for different values 
of mff c , ct , and pr- The efficiency, which is tabulated 
in Ref. [25|, reaches a maximal value of e = 52% for 
A —> n + n~ decays with to = 2GeV/c 2 , pt > 4GeV/c, 
and ct = 6 cm. The dominant factor affecting e is the 
average transverse flight distance (r) = ct ( pr) /to. Re¬ 
flecting the 1 < r < 50 cm requirement, e drops rapidly 
when (r) goes below 1 cm or above 50 cm. In addition, 
e has some dependence on the A polar-angle 9 , mea¬ 
sured with respect to the direction of the e + e _ center 
of mass. For a 1 + cos 2 9 distribution in the CM frame, 
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the strongest dependence is observed for track momen¬ 
tum p < 0.3 GeV/c, where e is decreased by 22% relative 
to that of a uniform cos 6 distribution. For p > 2 GeV/c, 
e varies by no more than 8%. Similarly, the efficiency 
depends weakly on whether L is a scalar or a vector par¬ 
ticle. For a longitudinally polarized vector, e typically 
varies by a few percent relative to the scalar case, with 
the greatest impact being an efficiency reduction of 25% 
for pt < 0.3 GeV/c, m = 7 GeV/c 2 . 

The dominant source of background consists of 
hadronic events with high track multiplicity, where large- 
do tracks originate mostly from Ji°, A, K ± 1 and de¬ 
cays, as well as particle interactions with detector mate¬ 
rial. Random overlaps of such tracks comprise the ma¬ 
jority of the background candidates. 

We extract the signal yield for each final state as a 
function of L mass with unbinned extended maximum- 
likelihood fits of the to distribution. The procedure is 
based on the fact that signal MC events peak in m while 
the background distribution varies slowly. The fit proba¬ 
bility density functions (PDFs) for signal and background 
are constructed separately for each mode and each data 
sample. The PDFs account for the signal mass resolu¬ 
tion, which is evaluated separately in each of 11 mass re¬ 
gions, where each region straddles the mff c value of one 
of the signal MC samples of the first type. In region i, the 
value of the signal PDF for a candidate with hypothesis 
mass mo, measured mass m, and mass resolution uncer¬ 
tainty a m is Pg(m) = H l s ((to — mo)/cr m ), where H l s (x ) 
is the histogram of the mass pull x = (m MC — m,^ 10 ) /a^ c 
for signal MC events of true mass m^ LC , measured mass 
m MC , and m MC uncertainty a^ G . This PDF accounts 
correctly for the large variation in a m with r and to. 

The background PDF Pb(to) is obtained from the 
data, so as not to rely on the background simulation, 
with the following procedure. First, we create a variable- 
bin-width histogram Hjjljn) of the data m distribution. 
The width of a histogram bin, whose lower edge is in 
to region i, is uy = nRi , where n = 15, and Ri is the 
RMS width of the signal to — TOq 10 distribution in that 
region. The value of Ri ranges from about 6MeV/c 2 for 
to“ c = 0.5 GeV/c 2 to 180 MeV/c 2 for toq 10 = 9.5 GeV/c 2 . 
We obtain Pb(to) by fitting I7 d(to) with a second-order 
polynomial spline, with knots located at the bin bound¬ 
aries. Simulation studies of the background mass distri¬ 
bution show that the choice n = 15 is sufficiently large 
to prevent Pb(to) from conforming to signal peaks and 
thus hiding statistically significant signals, yet sufficiently 
small to avoid high false-signal yields due to background 
fluctuations. Fig.[T]shows the to distributions of the data 
(with uniform mass bins) and the background PDFs. 

We scan the data in search of an L signal, varying 
too in steps of 2 MeV/c 2 . At each scan point, we fit the 
data in the full mass range using the PDF ngPg + tibPb, 
where the signal and background yields ng and ns are 
determined in the fit. The statistical significance S = 



m [GeV/c 2 ] 

FIG. 1: Mass distribution of the T(4S) + off-resonance data 
(red solid points) and Y(3S) + T(2S) data (blue open squares) 
for each mode, overlaid with the background PDF Pb in 
matching color. In the fi + mode, the bin width in the 
range m < 370MeV/c 2 is 10MeV/c 2 . 


sign (n s )y / 2log(£s7^B), where Cg is the maximum like¬ 
lihood for n s signal events over the background yield, 
and Cg is the likelihood for ng = 0, is calculated for 
each scan point. The distributions of S values for all the 
scan points are nearly normal. 

Significance values greater than 3 are found in two 
scan points, both in the /r + /i~ mode in the T(4S) + 
off-resonance sample. The highest significance is S = 4.7, 
with a signal yield of 13 events at the low-mass threshold 
of too = 0.212 GeV/c 2 . The second-highest significance 
of S = 4.2 occurs at mo = 1.24 GeV/c 2 , corresponding 
to a signal yield of 10 events. To obtain the p -values 
for these significances, we perform the scans on a large 
number of spectra generated according to the 

background PDF, obtained from the data with a finer 
binning of n = 5. With this choice of n, the generated 
spectra are not sensitive to fluctuations of the order of 
the signal resolution (which correspond to n = 1), yet 
include features that are much smaller than the resolu¬ 
tion of the PDF (n = 15). We find that the probability 
for S > 4.7 (4.2) anywhere in the spectrum with 

< 0-37 GeV/c 2 (to p + m - > 0.5 GeV/c 2 ) is 4 x 10 -4 
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(8 x 10" 3 ). The p-values are consistent with the naive 
expectation p(S)w/R , where p(S) is the p- value without 
the “look-elsewhere effect”, w is the width of the mass 
region under study, and R is the average value of R t . We 
do not include the other modes in the calculation of the 
p-values. Doing so would naively multiply the p-values 
by about six. Further study provides strong indication 
for material-interaction background in the 0.212 GeV/c 2 
region. Specifically, most of the 34 p, + p,~ vertices with 
m< 0.215 GeV/c 2 occur inside or at the edge of 
detector-material regions, including 10 of the vertices 
that also pass the e + e _ selection criteria and 10 that 
pass the criteria. Thus, the peak is consistent 

with misidentified photon conversions and hadronic in¬ 
teractions close to the mass threshold. We conclude that 
a significant signal is not observed. 

Systematic uncertainties on the signal yields are cal¬ 
culated for each scan fit separately. The dominant un¬ 
certainty is due to the background PDF, and is evalu¬ 
ated by repeating the scans with n = 20, which is the 
maximal plausible value for n that does not lead to a 
large probability for false-signal detection. This uncer¬ 
tainty is a few signal events on average, and generally 
decreases with mass. An additional uncertainty is eval¬ 
uated by taking f?,; from events with pt < 0.8 GeV/c or 
Pt > 0.8 GeV/c. To estimate uncertainties due to the 
weak signal PDF dependence on r and m, we repeat the 
scans after obtaining H l s from signal MC events with ei¬ 
ther r < 4 cm or r > 4 cm, as well as from signal MC 
events from adjacent mass regions. The uncertainty due 
to the signal mass resolution is evaluated by comparing 
the mass pull distributions of K® mesons in data and 
MC, whose widths differ by 5%. A conservative uncer¬ 
tainty of 2% on the signal reconstruction efficiencies is 
estimated from the K ° reconstruction efficiency in data 
and MC. Smaller uncertainties on the efficiency, of up 
to 1%, arise from particle identification, and signal MC 
statistics. The total uncertainties on the efficiency are 
reported in the efficiency tables i25]. 

Observing that the likelihood Cs is a nearly normal 
function of the signal yield, it is analytically convolved 
with a Gaussian representing the systematic uncertain¬ 
ties in ns, obtaining the modified likelihood function C' s . 
The 90% confidence level upper limit Us on the signal 
yield is calculated from f^ s C'dns // 0 °° C'dns = 0.9. Di¬ 
viding Us by the luminosity yields an upper limit on the 
product er(e + e _ -4 LX) B(L —»• f) e(f). This limit is 
shown for each mode as a function of mo in Fig. [2] and 
given in the supplemental material [25| . 

Determining the efficiency from the B —>• X S L signal 
MC sample, we obtain upper limits on the product of 
branching fractions B(B —» X S L)B{L —> f) for each of 
the final states /. These limits are shown in Fig. [3] 

In conclusion, we have performed a search for long- 
lived particles L produced in e + e _ collisions. No signal 
is observed, and upper limits on cr(e + e _ —> LX) B(L —> 



2X10' 1 1 2 345678 10 

m 0 [GeV/c 2 ] 


FIG. 2: The 90% confidence level upper limits on <r(e + e _ —► 
LX)B(L —> f)t(f) as a function of L mass for the T(45) + 
off-resonance sample (red lower points) and for the T(3S) + 
T(2S) sample (blue upper points). The limits include the 
systematic uncertainties on the signal yield. 

/) e(/) and on B(B -A X S L)B{L —> f) are set at 90% 
confidence level for six two-body final states /. We pro¬ 
vide detailed efficiency tables to enable application of our 
results to any specific model j2f|. 
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FIG. 3: Implications of the results for Higgs-portal scenarios, 
showing the 90% confidence level upper limits on the product 
of branching fractions B(B —> X a L)B{L —> /) as a function 
of L mass for each final state / and for different values of ct. 
The limits include all systematic uncertainties. 
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